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ABSTRACT 

A family of conditionally stable, forward Euler finite-difference equations can be constructed for 
the simplest equation of Schrodinger type, namely u, = iu«. Generalization of this result to physically 
realistic Schrodinger type equations is presented. 
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CONSTRUCTION OF STABLE EXPLICIT FINITE-DIFFERENCE 
SCHEMES FOR SCHRODINGER TYPE DIFFERENTIAL EQUATIONS 


Introduction 


Discrete finite-difference models of differential equations have been a traditional and popular 
technique for obtaining numerical solutions of both ordinary and partial differential equations [1,2]. 
An essential component of this procedure is the replacement of the derivative by its discrete analogue 
based on the definition of the derivative as a limit process; that is 

dx T . x(t +h) - x(t) 

ar ■ J — fc — . (i) 

n-*o 

and 


dx x* +1 - x* 
It H ’ 


( 2 ) 


where x* is an approximation to x(U) with Ik = hk. However, more general definitions of the derivative 
can be defined. As an example consider 


dx 

dt 


= Lim 

h-0 


X(t +h) ~ X(t) 

m 


(3) 


where 0(h) has the property 


0(h) = h +0(h 2 ). 

Note that for this case the discrete replacement for the derivative is 


(4) 


dx x*+i - x* 

(Jt 0(E) ’ 


(5) 


which, in general, differs from Eq. (2) for finite values of the step-size h. 

The major implication of the above remarks is that new discrete models of differential equations 
can be constructed based on Eq. (3). These models have the potential for providing solutions that arise 
when the conventional discrete derivation of Eq. (2) is used. 

As an illustration, consider the differential equation 

= rtx, x(to) = xo, A = real parameter, (6) 


dx 

It 


with solution 


48 


f 


( 7 ) 


x(t) = xo exp [ U (t — to) ]. 

Application of conventional procedures [1] gives, for example, the finite-difference scheme 

-Bi = iAx*. ( 8 ) 

n 

In contrast, consider the following exact scheme [3] 

— — = Ux*, xo = given, ( 9 ) 

where 


m = 



( 10 ) 


An easy and direct calculation [4] shows that for Eqs. (6) and (8), we have for their solutions 
x(U) * xfc, while for Eqs. (6) and (9), x(t*) = x* for all step-sizes h. Also, from Eq. (10), we obtain 


0(h) = h + iAh 2 + 0(h 3 ). 


( 11 ) 


Consider the simplest Schrodinger-type equation 


Conventional finite-differences, using a forward Euler replacement for the time-derivative, give a 
scheme that is unconditionally unstable [5J. We now show how conditionally stable schemes can be 
constructed. 


To proceed, based on the discussion presented earlier, we construct the following finite-difference 
scheme 


Um +1 — U^i _ Um +/ 2ll/n + Um -7 (13) 

101 ( 50 ) 02 ( 50 ) 

where the “denominator functions” have the properties 

0i(At, A) =At + iA(At) 2 + 0[(At) 3 ] (l 4a ) 

02 (Ax,t) = (Ax) 2 + 0[(Ax) 3 ] (l 4b ) 

and A, for the moment, is an unspecified parameter. Defining R(At, Ax, A) as 

R = =Ri(At,Ax,A) + iR 2 (At, Ax, A) (15) 

we can rewrite Eq. (13) in the form 
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( 17 ) 


i4 +i = RuS, + i + (1 -2R) v&n + Ru£. 

The substitution of a typical Fourier mode 

uH, = C(n) exp [ ky(Ax)m], 

into Eq. (16) and requiring the C(n) be bounded allows the determination of the stability properties. 
(This concept of stability is called “practical stability” [5,6].) A straight-forward calculation shows that 
Eq. (16) is stable if the following condition is satisfied. 

( R i - + ( r 2) 2 * jg • (18) 

This relation has the interesting geometric interpretation: In the (Ri,R 2 ) plane, Eq. (16) is stable for 
points on and inside the circle of radius 0.25 centered at (0.25, 0). We will refer to the inequality of Eq. 
(18) as the circle condition. 

For a given application, the following procedure is to be followed: 

a) Select denominator functions with the properties given in Eqs. (14) and calculate R in Eq. (15). 

b) Choose a point (Ri, R 2 ) consistent with the circle condition of Eq. (18). 

_ c ) Determine Ri(At, Ax, A) and R 2 (At, Ax, A) from Eq. (15) and set them equal, respectively to 
Ri and R 2 , i.e., 

Ri(At, Ax, A) =Ri, R 2 (At, Ax, A) = R 2 . (19) 

d) Select a value for the space-step, Ax, and solve the two relations of Eq. (19) for At and A in 
terms of Ax. Doing this gives 

At = fi(Ax), A =f 2 (Ax). (20) 

Thus, the selection of the point (Rj, R 2 ) , satisfying the circle condition and the relations of Eq. (20), 
defines completely the finite-difference scheme of Eq. (16). 

As an example, consider the following denominator functions 

<t>\ = At + iA(At) 2 , <p 2 = (Ax) 2 , (21) 


with 


R = 


1 + i 
4 ' 


The following results are easily obtained 


( 22 ) 
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(23) 


A, = ^, A = — 

and 

* +1 = + >*-1] + <24) 

Note that everything can be expressed in terms of Ax, the space step-size; its value can be selected as 
desired. 

A second, nontrivial example is the nonlinear Schrodinger equation 

U< = iUjcr + | u 1 2 u. (25) 

A finite-difference scheme that embodies the work of this paper and also uses the nonlocal modeling 
of the nonlinear term [3] is 

Um + l ~ 2u ot + Um-1 

<P 2 


Um +1 


- Um 


<(>i 


= i 



lCi + 1 + Um-1 

♦ 

Um + 1 

+ Um-1 

2 



2 


„« + l 
Um > 


(26) 


where the star (*) denotes complex conjugation, and the denominator functions ^i(At, A.) and 
0 2 ( Ax, A) satisfy the conditions of Eqs. (14). Note that Eq. (20) is an explicit finite-difference scheme. 


In summary, we have presented a new procedure to construct explicit finite-difference schemes 
for Schrodinger type partial differential equations. In general, we expect these schemes to be condi- 
tionally stable. This paper summarizes research that has already been published in Physics Review A , 
June 1989 and in the “Proceedings of the 2nd IMAC Conference on Computational Acoustics.” 
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